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Prepotential approach to exact and quasi-exact solvabilities of 
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In this talk I present a simple and unified approach to both exact and quasi-exact 
solvabilities of the one-dimensional Schrodinger equation. It is based on the prepotential 
together with Bethe ansatz equations. This approach gives the potential as well as the 
eigenfunctions and eigenvalues simultaneously. In this approach the system is completely 
defined by the choice of the change of variables, and the so-called zero-th order prepo- 
tential. We illustrate the approach by several examples of Hermitian and non-Hermitian 
Hamiltonians with real energies. The method can be easily extended to the constructions 
of exactly and quasi-exactly solvable Dirac, Pauli, and Fokker-Planck equations, and to 
quasinormal modes. 
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1. Introduction 

I am greatly honored to be able to participate in this conference celebrating Prof. 
Yang's 85th birthday. We wish Prof. Yang a very happy birthday, and many more 
to come. 

In this talk I would like to present a constructive approach to both exact and 
quasi-exact solvable one-dimensional Schrodinger equations. Everyone knows what 
exact solvability means. It is fair to say that all of us learned the principles of 
quantum mechanics through several well-known exactly solvable (ES) models, such 
as the infinite square well, the harmonic oscillator, the hydrogen atom, etc. But 
in actual life, exactly solvable systems are hard to come by. Most systems one 
encounters are non-solvable, and one has to adopt various approximation schemes 
to solve them. 

Two decades ago, a new class of potentials which are intermediate to ES poten- 
tials and non solvable ones have been found for the Schrodinger equation. These 
are called quasi-exactly solvable (QES) models for which it is possible to determine 
algebraically a part of the spectrum but not the whole spectrum i>2.3,4,5,6 r^j^^ 
plest QES model first discovered is the sextic oscillator ^ . The discovery of this class 
of spectral problems has greatly enlarged the number of physical systems which we 
can study analytically. 

Here I would like to discuss a novel approach that treats both ES and QES 
systems on the same footing ^. It is based on the so-called prepotential together 
with the Bethe ansatz equations. 
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Since QES systems are less familiar to most people, I shall give a very brief 
introduction to QES theory before the prepotential approach is discussed. 



2. A brief introduction to QES theory 

The essence of quasi-exact solvability is most easily explained as follows ^. It is 
well known that any Hamiltonian H can be represented as an infinite-dimensional 
Hermitian matrix 



(1) 
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where the elements Hnm = (V'nl^^lV'm) depend on the choice of a complete set of 
orthonormal functions (,i'n forming a basis of the Hilbert space. The solution of the 
spectral problem then reduces to a diagonalization of matrix {-ff„,m}- Unfortunately, 
contrary to the case of finite matrices, there is no general algebraic rules that would 
allow one to diagonalize the infinite-dimensional H. 

If the matrix {Hnm} is very specific so that it can be reduced to the diagonal 
form with the aid of an algebraic process, then the quantal system is exactly solvable. 



H 



HoQ 

Fii 

^22 

H33 



(2) 



The most well-known example is the harmonic oscillator. 
Suppose a Hamiltonian H is reducible to a block form 
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(3) 



where the block in the upper left corner is an {N -|- 1) x {N + 1) Hermitian matrix, 
and the block with asterisks is an infinite matrix with non-vanishing elements. The 
upper left block can be diagonalized without touching the infinite one. This means 
that one can determine only a part of the spectrum of H with A'' -|- 1 eigenvalues 
and eigenfunctions. This system is called QES. 
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QES models are not simply mathematical constructs. Although the first and 
simplest example of QES model is the sextic oscillator, physical QES systems have 
been found. For examples, the system of two electrons in a three-dimensional oscilla- 
tor potential ^, two-dimensional charged particle in Coulomb and uniform magnetic 
fields (the Schrodinger, Klein-Gordon, and Dirac cases are discussed in Refs. 9, 10 
and 11, respectively), and electron on a lattice in magnetic field (Azbel-Hofstadter 
problem) 

Usually a QES problem admits a certain underlying Lie algebraic symmetry 
which is responsible for the quasi-exact solutions. Such underlying symmetry is 
most easily studied in the Lie-algebraic approach The essence of this approach 
is as follows. Consider a Schrodinger equation Htjj = Eij} with Hamiltonian H = 
—(P/dx^ + V{x) and wave function 4'{x). Here x belongs either to the interval 
(— oo, oo) or [0, oo). Now suppose we make an "imaginary gauge transformation" on 
the function (j): (j){x) = (f>{x)e~^°^^\ For physical systems which we are interested 
in, the phase factor exp(— W^o(a;)) is responsible for the asymptotic behaviors of 
the wave function so as to ensure normalizability. The function ^(a;) satisfies a 
Schrodinger equation with a gauge transformed Hamiltonian Hw = e^°He~^° . 
Suppose Hw can be written as a quadratic combination of the generators J" of 
some Lie algebra with a finite dimensional representation, i.e., 

Hw = ^ CabJ^J^ + ^ Co J" real constant, (4) 

where the constants Cah and are usually taken to be real 3.5,6 q^]^gj^ within 
this finite dimensional Hilbert space the Hamiltonian Hw can be diagonalized, and 
therefore a finite number of eigenstates are solvable. Hence the system described 
by H is QES. For one-dimensional QES systems the most general Lie algebra is 
sl{2). The generators J° of the sl{2) Lie algebra take the differential forms: J+ = 
z^dz — nz , = zdz — n/2 , J~ = dz [n = 0, 1,2,.. .). The variables x and z 
are related by some function, n is the degree of the cigcnfunctions (i, which are 
polynomials in a (n + l)-dimensional Hilbert space with the basis (1, z, z^, . . . , z"). 

The Lie-algebraic approach to QES models excels in revealing the underlying 
symmetry of a QES system explicitly. However, solutions of QES states are more 
directly found in the analytic approach based on the Bcthe ansatz equations ^. In 
this analytic approach the form of the wave functions containing some parameters 
are assumed from the very beginning, and these parameters are fitted to make the 
ansatz compatible with the potential under consideration . 

Recently, a different approach to QES models has emerged i3,i4 ^Yas approach 
the emphasis shifts from the potential to the so-called prepotential (or superpoten- 
tial), a concept which plays a fundamental role in supersymmetric quantum me- 
chanics (wc note here that the prepotential considered here is the integral of the 
superpotetial in Refs. 13 and 15). Prepotential has been extensively employed to 
study classical and quantum integrability in Calogero-Moser systems The merit 
of this approach is that the form of the potential of the system concerned needs 
not be assumed from the beginning. All information about the system is contained 
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in the prepotential and the solutions, or roots, of the Bethe ansatz equations (in 
terms of certain transformed coordinates). The prepotential and the roots determine 
the potential as well as the eigenfunctions and eigenvalues simultaneously. Also, in 
this approach exact and quasi-exact solvabilities can be treated on the same foot- 
ing. Furthermore, such approach facilitates extension of the QES theory from the 
Schrodinger equation to equations for multi-component wave functions 1^,23 More 
recently, QES theory was extended to the Fokker-Planck equations also via the 
prepotential approach 

The emphasis of the works just mentioned was placed on the feasibility and 
elegance of the prepotential approach. But in these works the forms of the prepo- 
tential and the required change of coordinates were either directly adapted from the 
known ES and QES models, or given as known for the new QES systems. To make 
the prepotential approach a satisfactory approach one must be able to determine the 
choice of coordinate transformation and the prepotential from the very beginning. 
This was recently achieved for certain classes of coordinates and prc^potentials ^. 
In what follows we shall present the main ideas of this approach, and demonstrate 
how the simplest ES and QES models could be easily derived. 

3. Prepotential approach 

Suppose (f)o{x) (— oo < X < oo)) is the ground state, with zero energy, of a Hamilto- 
nian ffo: ^?o0o = 0. By the well known oscillation theorem cpo is nodeless, and thus 
can be written as (po = e^'^oix) ^ where Wo(a;) is a regular function of x. For the 
squarc-integrable (f>Q, this is the simplest example of quasi-exact solvability. This 
implies that the potential Vq is completely determined by Wq: Vq = Wq — Wq , and 
consequently, the Hamiltonian is factorizable (we adopt the unit system in which h 
and the mass m of the particle are such that H = 2m =1): 

We shall call Wq{x) the zero-tli order prepotential. 

Consider now a wave function {N > 0) which is related to of Ho by 
<i>N = (po4>N, where 

4)N = {z- zi){z - Z2)- ■ -{z - Zn), ^0 = 1- (6) 

Here z = z{x) is some function of x. In taking the form of (1>n in Eq. (6) I have 
assumed that the only singularities of the system are z = ±oo. Other situations 
are discussed in Ref. 7. The function ^jv is a polynomial in an [N + l)-dimensional 
Hilbert space with the basis {l,z,z'^,. . . ,z^). One can rewrite (j)N as 

= exp {-Wn{x, {zk})) , (7) 

with the N-th order prepotential Wn being defined by 

N 

WN{x,{zk}) = Wo{x) -J2^u\z{x) - Zk\. (8) 
fe=i 
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Operating on (j)N by the operator —(P/dx^ results in a Schrodinger equation 
Hn4'n = 0, where 

Vn = W'j^ - W^. (10) 
From Eq. (8), Vn has the form Vn = Vo + AVn, where Vq = W'^ - W^', and 

^v„.-2(«.i,--^)|:j2_+5^_|_. ,„) 

Here the prime denotes derivative w.r.t. the variable x. 

AVat is generally a meromorphie function of z with at most simple poles. Let 
us demand that the residues of the simple poles, Zk, k = 1,. . . ,N should all vanish. 
This will result in a set of algebraic equations which the parameters {zk} must 
satisfy. These equations are called the Bethe ansatz equations for {zk}- With {zk} 
satisfying the Bethe ansatz equations, AVn will have no simple poles at {zk} but 
it still generally depends on {zk}- Thus the form of Vn is determined by the choice 
of z'2 and Wo{x). 

In what follows we would like to demonstrate that the choice of z'^ and Wqz' 
determine the nature of solvability of the quantal system. We shall restrict our 
consideration only to those cases where Wqz' = P„i{z), z'^ = Q„(z) and z" ~ Ri(z) 
are polynomials in z of degree m, n and I, respectively. Of course, Qn{z) and Ri{z) 
are not independent. In fact, from 2z" = dz''^/dz we have I = n — 1 and 

z _i?„_i(^)_-^_. (12) 

Consequently, the variables x and z are related by (we assume z{x) is invertible for 
practical purposes) 

and the prepotential Wo{x) is determined as 



Wo{x) = J 




(14) 

z=z{x) 



Eqs. (13) and (14) define the transformation z{x) and the corresponding prepo- 
tential Wo(x). Thus, Pm{z) and Qn{z) determine the quantum system. Of course, 
the choice of Pm and Qn must be such that Wq derived from (14) must ensure 
normalizability of 0o = exp(— Wo). 

Now depending on the degrees of the polynomials Pm and Qn, we have the 
following situations: 
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(i) if max{m, n — 1} < 1, then in Vn{x) the parameter N and the roots Zk's wiU 
only appear as an additive constant and not in any term involving powers of z. 
Such system is then exactly solvable; 

(ii) if max{m. n — 1} = 2, then N will appear in the first power term in z, but z^'s 
only in an additive term. This system then belongs to the so-called type 1 QES 
system defined in Ref. 3, i.e., for each A'' > 0, Vjv admits A'' + 1 solvable states 
with the eigenvalues being given by the N + 1 sets of roots Zk's. This is the 
main type of QES systems considered in the literature; 

(iii) if min{m, n — 1} > 3, then not only N but also Zk's will appear in terms 
involving powers of z. This means that for each N >0, there are A'' + 1 different 
potentials V^V; differing in several parameters in terms involving powers of z, 
have the same eigenvalue (when the additive constant, or the zero point, is 
appropriately adjusted). When z^s appear only in the first power term in z, 
such systems are called type 2 QES systems in Ref. 3. We see that QES models 
of higher types are possible. 

These general situations are considered for cases where z'^ = Q2(z) are quadratic 
in z (this includes Qi and Qo as special cases if some of the coefficients vanish) in 
Ref. 7. This choice of z'^ covers most of the known ES shape-invariant potentials 
in Ref. 15 and the s?(2)-based QES systems in Ref. 3, and a new one discussed in 
Refs. 14 and 18. Such coordinates are called "sinusoidal coordinates", which include 
quadratic polynomials, trigonometric, hyperbolic, and exponential types. The con- 
nection of sinusoidal coordinates with ES theory has been extensively discussed in 
Ref. 19. For simplicity, in this talk I shall only discuss cases in which z'^ is linear 
in z, i.e., z'^ = Qi(z) = qiz + qo, where qi and go are real constants. 

With this choice of z'^, we have 

In deriving Eq. (15) use has been made of the following identities: 

^ 1 1 / 1 \ 

£ - ^k){z - Zi) '^^^ Z- Zk\Zk- zi) ' ^^^"^ 

N ^ 1 / \ 

Et Vr V = 2E^f^V (17) 

ifj^i - Zk){z - Zl) Z-Zk \Zk-Zlj 

kjtl k^l 

Demanding the residues at Zfc's vanish gives the set of Bethe ansatz equations 

^-(^O-T-E^^=0' k = l,2,...,N. (18) 

Putting back the set of roots Zk into Eq. (15), we obtain a potential Vn{x) without 
simple poles. The degree of Pm{z) determines the nature of the solvability of the 
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system, namely, for m — 1,2,3,..., the system is, respectively, ES, type 1 QES, 
type 2 QES, and higher types QES, as discussed generally in the last section. 

To proceed further, we must specify Pm{z)- Below I shall discuss cases with 
m < 2, i.e., only ES and type 1 QES models. 



4. Examples 

Let us take Qi{z) = AAz + qo {qi = 4A), where A, go are real constants. This 
implies z"{x) = Ro{z) = 2A. Hence, the general solution of z{x) is a quadratic form 
of x: 

z{x) = Ax^ -\-Bx-\-C, A, B, C: real constants. (19) 

It is easily checked that qo is related to A, B and C through z'"^ = AAz + B"^ —AAC. 

Wc now show how some ES and QES models can be constructed in the prepo- 
tential approach by taking different values of m. 



4.1. Exactly solvable cases: m = 1 

Suppose m = 1 so that Pi{z) = Aiz + Aq {Ai, Aq real). By writing Pi = A-i{z 
Zk) + AiZk + ^0, one obtains from Eq. (15) that 

Ay;v = -2AiVl-2V^|A(z,)-^-E^^ 

= -2AiN. (20) 

The last term in braces in Eq. (20) vanishes when ^^'s satisfy the Bethe ansatz 
equations (18). Now N only appears as a parameter in an additive term, and not 
in terms involving powers of z in AV/v- The roots z^s do not appear at all. The 
additive term can be treated as the eigenvalue. The Schrodinger equation reads 

[--^+ W'^ - W'^ e-^~ = 2AiiVe-^«. (21) 

We see that the potential Wq — Wq is ES: by varying N ^ one obtains all the 
eigenvalues 2A-yN and the eigenfunctions = exp(— VFat). 

As an example, let us take Q\{z) = 1 and P\{z) = bz. A particular solution of 
z is z{x) = X. Prom Eq. (14) one gets 

r hx bx"^ 

Wq{x) = / —dx = — + const. (22) 

The integration constant is to be determined by normalization of the wave function. 
For (f)Q = exp(— Wo) to be square-integrable, one must assume b > 0. The BAE (18) 
are: 

hxk-Y — ^- — =0, k = l,...,N, (23) 
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and the potential is 

Vo = b^x^ - b, AVn = -2Nb. (24) 
This leads to the Schrodinger equation: 

(^-^+ b^x^^ e-^« = b{2N + l)e-^". (25) 

This system is just the well-known simple harmonic oscillator. 

Wc note here that by rcscaling Vbxk Xk, Eq. (23) will have 6 = 1. The 
resulted equations are the equations that determine the zeros of the Hermite poly- 
nomials Hn{x) as found by Stieltjes 20,21,16 jjguQg .^^g have reproduced the well 
known wave functions for the harmonic oscillator, namely, (pN = exp(— VFjv) ~ 
ex.p{-bx^/2)HNiVbx). 

4.2. Type 1 quasi- exactly solvable cases: m — 2 

Next we consider P2{z) = ^2^^ + Aiz + Aq. By a similar argument we obtain 

I^Vn = -2A2Nz-2A2^Zk-2AiN. (26) 
fe=i 

The Schrodinger can be written as 

(-^ + W;,^ - < - 2A2N^ e-^^ = 2 (^^2 5^ 2fe + A^N^ e"^". (27) 

Unlike the previous case, now N not only appears in an additive constant term but 
also in the term with z, and the set of roots z^'s appear in the additive term. This 
system is the so-called type 1 QES models. Type 1 QES models classified as class 
VI in Ref. 3 belong to this category. 

A well-known example is the sextic oscillator, the simplest QES model of this 
type ^, which is termed as class VI in Ref. 3. In our prepotential approach, this 
system is defined by z{x) = .x^ and ^2(2) — 2{az^ + bz). Then Qi{z) = Az, and 

r ax^ + bx^ , 1 4 1,0 /„ox 
Wq{x) = / = — dx = -ax^ + -bx^ + const. (28) 



4 2 

Here a > to ensure square-integrability of the wave function. The BAE (18) are: 
2azl + 2bzk - I - 4:y^ =0, k = l,...,N, (29) 

and the potential is 

Vn = a^x^ 2abx'^ ^ [6^ - (47V ^ 3) a] x'^-AoY^z^- {AN + 1)6. (30) 

k 

This is just the well-known QES sextic oscillator first discovered in Ref. 1. 

These two examples show how simply and directly the two most classic repre- 
sentatives of ES and QES models are derived in the prepotential approach. Other 



2, 2008 16:38 WSPC/INSTRUCTION FILE HoCL-CNYang85 



9 



ES and QES models can be constructed in a similar way. Particularly, we can con- 
sider cases where z'^(x) = Q2{z) with q2 ^ 0. Here z{x) is again some sinusoidal 
coordinates, which include the exponential, hyperbolic and trigonometric functions. 
By taking appropriate P„i{z), one can reconstruct class I, II and X QES models in 
Ref. 3, and some of the ES models listed in Ref. 15, namely, the Morse potential, 
the Scarf I and II potentials, and the Poschl- Teller potential ^. 



5. Non-Hermitian potentials with real spectra 

Now we turn to a different kind of spectral problem. About ten years ago, it was 
noted that with properly defined boundary conditions the spectrum of the non- 
Hermitian Hamiltonian H = + x'^{ixY {u > 0) is real and positive It is 
believed that the reality of the spectrum of this Hamiltonian is a consequence of 
its unbroken VT symmetry. This discovery immediately sparked great interest in 
searching for new systems with such properties, and in understanding the reasons 
behind the reality of their eigenvalues (for review, see e.g. Ref. 23). 

Soon after the appearance of Ref. 22, it was found that, by allowing non- 
Hermitian "PT-symmetric Hamiltonians, a QES polynomial potential can be quartic 
in its variable Previously, it was believed the lowest-degree one-dimensional QES 
polynomial potential is sextic, namely, that presented in Sect. 4. Now let us demon- 
strate how simply the Hamiltonian discovered in Ref. 24 is derived in our approach 
if we allow Pm{z), or equivalently Wo{x), to be complex. 



5.1. VT-symmetric quartic case 

We take z = x (i.e. Qo = 1), and P2{x) = ax^ + hx -\- c for type 1 QES system. 
Eq. (14) then gives Wq{x) = ax'^/3 + bx'^/2 + ex. It is obvious that, with a, b and 
c real, the wave functions are not normalizable as the factor exp{—Wo) does not 
converge on the whole line. Hence a QES Hamiltonian with quartic polynomial 
potential is impossible with real coefficients in Wq. 

What if we allow P2{x) to be complex? Suppose we take a = ia, b = (3 and 
c = t7, where a, /? > and 7 are real. Obviously exp(— Wq) is now square-integrable 
on the whole line for 6 > 0. Prom Eqs. (15) and (18) we obtain straightforwardly 
the type 1 QES potential 



Vn{x) = -a^x'^ + 2ial3x^ -\- {0^ - 2a-i) x^ - 2i [{N + l)a - (3j] x 
^'^ + l3 + 2ia^Xk 



with the BAE 



iaxl + I5xk + i"i — ^ — = 0, fc = 1, 2, . . . , iV. 

^ Xk-Xl 



l^k 



(31) 



(32) 



Vn given by Eq. (31) is precisely the 'PT-symmetric QES quartic potential obtained 
in Ref. 24 (in the notation of Ref. 24, the potential has parameters a = 1, /3 = 
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a, J = b, and N + 1 = J). In Ref. 24 it was shown that this T-'T-symmetric 
Hamiltonian is also related to the Lie-algebra sl{2, C). It is nice to see that in our 
prepotential approach this potential can be so simply and directly derived, without 
any knowledge of its underlying symmetries. 

The properties of this system have been fully described in Ref. 24, so we will 
not attempt a full discussion of it. Let us just consider the two simplest cases with 
= and = 1. 

For A'' = 0, the potential is just Vq with only one solvable state corresponding 
to energy E = Q. For N = 1, there are two solvable states. Prom Eq. (31) we can 
take the potential to be 

Vi = -o?x^ + 2ia^x^ + (/J^ - 2a7) - 2i{2a - /?7)x (33) 

with eigenvalues 7^ + /? + 2iax\. The roots xi satisfy the BAE (32): iax\ + fixi + 
Z7 = 0. The solutions are: xi = i{(3 ± ^ 0'^ + Aa^)/2a. Hence we see that for 
/32 + 4a7 > the non-Hermitian potential Vi has two normalizable states with real 
energies 7^ ± \/ (3"^ + 'iaj. But if + 4a7 < 0, the energies are complex. This is 
consistent with the numerical results in Ref. 24. For higher values of N the analysis 
is more complicated and numerics are needed 

5.2. Sextic case without T"T symmetry 

We now consider a second non-Hermitian QES system. Suppose we assume the 
parameter b in the sextic oscillator (class VI in Ref. 3) in Sect. 4.2 to be purely 
imaginary: b = ip, where (3 and a > are real. The factor exp(— Wo) is square- 
integrable. Now the potential and BAE are: 

Vat = a^x^ + 2ial3x^ - [0^ + {AN + 3) a] x^ - Aa^Zk - i{AN + (34) 

k 

2azfe + 2i/3zfc - 1 -4 V = 0, k = l,...,N. (35) 

Again, this non-Hermitian system could have real eigenvalues. Trivially iov N = Q 

the potential Vo admits one normalizable cigenstate with real energy E = Q. For 
N = 2, the solutions of Eq. (35) are: zi = (-i/3 ± s/2a- j3'^)/2a. Then Eq. (34) 
becomes 

Vi = a^x^ + 2iaf3x^ - [0^ + la)x^ - 3z/3 ± 2s/2a- P"^. (36) 

Hence the non-Hermitian potential defined by the first four terms in Eq. (36) has 
two real eigenvalues E = ±2i/2a — /J^ if < 2a. The energies become complex if 
> 2a. Like the quartic case just discussed, for higher values of N the analysis of 
its spectrum requires the aid of numerics. Let us emphasize that this systems is not 
■PT-symmetric. 

We note here that the same analysis can be applied to class I and VII QES 
systems in Ref. 3 to get other non-Hermitian QES systems with real energies. 
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6. Summary 

I hope I had given you a general idea of how simple and elegant the prepotential 
approach to ES and QES models is. In this talk I have limited myself to discussing 
our approach only to systems defined on the full line owing to time limitation. Our 
approach has been extended to ES and QES systems defined on half-line and on 
finite intervals ^. It can be easily extended to the construction of exactly and quasi- 
exactly solvable Dirac, Pauli Fokker-Planck equations and to quasinormal 
modes, which are eigen-modes with complex energies of a Hermitian Hamiltonian 

25 

Of course, what I have presented here is quite preliminary. Our approach could 
be further developed. Firstly, so far we have only considered cases for which Wqz' 

and z'^ are polynomials in z. It is interesting to consider other possibilities for these 
two defining functions. Secondly, it is extremely interesting to see how our for- 
malism could be extended to many-body systems, such as the renowned Calogero- 
Sutherland-Moser systems. The analogues of Bethe-ansatz equations for these sys- 
tems are not known yet. Lastly, although some non-Hermitian Hamiltonians can 
be generated in this approach by making some coefficients in Pm{z) complex, as 
we have shown here, a full treatment of these kind of systems usually requires one 
to extend the basic variable x to the complex plane 22,23,24_ jj^jj^g ^ better under- 
standing of the prepotential approach to non-Hermitian Hamiltonians is needed. 
We hope to report on our progress in these directions in the near future. 
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